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This is a review paper for PAB3D’s history in the implementation of turbulence models 
for simulating jet and nozzle flows. We describe different turbulence models used in the 
simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations 
use modified linear or nonlinear two-equation models to account for supersonic flow as well 
as high temperature mixing. Two multiscale-type turbulence models are used for unsteady 
flow simulations. These models require modifications to the Reynolds Averaged Navier- 
Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two- 
equation (k-g) model with a RANS/LES transition function, dependent on grid spacing and 
the computed turbulence length scale. The second scheme is a modified version of the 
partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented 
in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational 
methods, code implementation, computed results for a wide range of nozzle configurations at 
various operating conditions, and comparisons with available experimental data. Very good 
agreement is shown between the numerical solutions and available experimental data over a 
wide range of operating conditions. 


I. Introduction 

K NOWLEDGE of jet mixing and nozzle aerodynamics is vital to several areas of commercial and military 
aircraft design, such as propulsion efficiency, propulsion integration, aeroacoustics, and jet interaction with 
aircraft structures. Initial jet flow conditions are determined by nozzle exit pressure, temperature, Mach number, and 
nozzle geometry. Once the exhaust flow leaves the nozzle, the jet flow forms a free shear layer. The action of 
turbulence dominates flow development farther downstream. As such, jet flow properties are difficult to measure or 
predict analytically. When Ludwig Prandtl introduced his mixing length hypothesis for turbulent flows 100 years 
ago, a brief analysis of a fully mixed jet was given as an example. For axisymmetric and two-dimensional jet flows, 
early analyses of mixing behavior were based mainly on this mixing length hypothesis and one-dimensional 
momentum theory (1-4). Mean flow properties for axisymmetric jets derived from these analytical models 
compared well with experimental measurements of jets at low subsonic speeds. However, data from jet flow 
measurements in the high subsonic and supersonic speed ranges (ref. 5) indicated significant departure from the 
results obtained by using one dimensional momentum theory. 

Jet and nozzle flow contain rich combinations of flow interactions and flow physics. These combinations include 
turbulent mixing as well as temperature and compressibility effects. Other factors may include chemical reactions or 
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shear layer instability. The main flow physics feature in the subsonic jet flow is shear layer development along the 
streamwise direction. The static pressure value is almost the same as the ambient pressure. In the absence of a 
pressure gradient, no significant inviscid flow feature will appear in a subsonic jet. According to reported 
experimental measurements, all turbulent axisymmetric subsonic cold jets below Mach 0.6 are similar if the flow 
variables are normalized by jet density and nozzle-exit velocity. 

On the other hand, supersonic jet flow features can be very complex. Because of the supersonic nozzle exit Mach 
number, a shock will form at the nozzle lip and the jet exit pressure can differ from the ambient pressure. This 
pressure difference between the jet and the ambient fluid must be resolved locally either across an oblique shock, by 
a prominent streamline curvature at the jet boundary, or both. A Mach disk can also form inside the jet. In addition, 
shocks initiated near the nozzle exit may reflect repeatedly at the sonic line in the shear layer downstream from the 
nozzle exit. Although the turbulence interacts with shocks in the jet, the position of the reflected shock depends 
mainly on the location of the sonic line in the turbulent shear layer. Such interdependence of flow interactions can 
become very complex. 

Earlier jet flow analysis codes, with or without chemical interactions included, were formulated with simplified 
assumptions of the Navier-Stokes equations and the turbulence model to provide the best jet flow simulation within 
modest limits of computing resources available during this time period. Analytical methods and simulation codes 
developed by this approach have been successfully applied to problems in air-breathing engine development, 
acoustics, and rocket propulsion. (6-12) However, there are some drawbacks to this approach. First, simplified 
assumptions are often difficult to justify. Second, application of the simplified formulations is limited to just the jet 
flow. These formulations are therefore difficult to integrate with computational codes for airframe aerodynamics 
when propulsion airframe integration analysis is required. It is preferable in such cases to perform the analysis with 
the three-dimensional Navier-Stokes equations without empirical assumptions for jet flow alone. 

Large temperature and pressure fluctuations have a profound effect on turbulence development in nozzle and jet 
flow. While several models have been developed to account for the effect of pressure fluctuations (compressibility 
correction models), very little has been done to account for large temperature fluctuations. This has led to poor CFD 
prediction of non-isothermal flows. For high- temperature jet flow, standard turbulence models lack the ability to 
predict the observed increase in growth rate of a mixing layer [13-14]. Several researchers [15-20] have modified 
one or more terms of the transport equations to obtain better agreement in high temperature flows. These 
modifications affect, directly or indirectly, the closure terms of the turbulent heat flux (pufi) and stresses 
( pUjU • ). Theis and Tam [15] changed several coefficients in the turbulent transport equations. However, such 

extensive modifications of model coefficients completely change the characteristics of the equations and may cause 
deficiencies in flow prediction accuracy, problems for which these turbulence models were originally designed. 
Other attempts to sensitize the turbulence model to temperature fluctuations involve more sophisticated closure for 
the turbulent heat flux term appearing in the average energy equation [16-18]. Explicit algebraic nonlinear heat flux 
models have also been tested for this purpose. These models have been successful in some fully-developed high 
temperature turbulent flows. A simpler approach was to model the value of as a function of the total temperature 
gradient in the flow [21]. 

For general use of steady- state jet flow simulation, some basic requirements must be met. The Navier-Stokes 
code should be upwind biased to capture shocks and other jet flow discontinuities. The code should also be fully 
three-dimensional in space because the relations between turbulent kinetic energy and Reynolds stresses are 
basically three-dimensional. The turbulence model should be capable of providing a time scale and a consistent 
description of the production and transport properties of turbulent kinetic energy. Therefore, at least a two-equation 
turbulence closure model is required. Many upwind-biased three-dimensional Navier-Stokes codes are available that 
meet the jet and nozzle flow simulation requirement. However, the availability of codes with a robust two-equation 
turbulence model in this class is limited. The PAB3D code has demonstrated the feasibility of establishing a unified 
method for subsonic and supersonic jet analysis with a general-purpose, three-dimensional Navier-Stokes code. 

PAB3D code was developed to obtain numerical solutions to the Reynolds Averaged Navier-Stokes equations in 
three-dimensional spatial domain. The main solver algorithm is the upwind Roe scheme, for which the numerical 
dissipation is small. The Jones-Launder (22) two-equation k - 8 turbulence closure model is used as the basis for all 
the corrections added to enhance simulation accuracy for a wide range of nozzle shapes and jet flows. This approach 
is chosen for jet flow analysis because it is consistent in tracking production and transport properties of turbulence 
kinetic energy and dissipation scale length in the shear flow. In the Jones-Launder k-z turbulence model, several 
empirical constants are required. Only the published values for these constants are implemented in the PAB3D code. 
These values are fixed for all computational applications of the PAB3D code. The final turbulence model 
formulation used in PAB3D contains the following major modifications: 
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1 - Near wall modification 

2- High temperature correction for hot jet 

3- Compressibility correction for supersonic mixing layer 

4- Algebraic Reynolds Stress for flow with separation 

This modified turbulence model reverts back to the standard k-c model in the flow regions where there is no need 
for any of these corrections. 

The limited capability of the Reynolds Averaged Navier- Stokes (RANS) approach, combined with eddy- 
viscosity turbulence models to simulate unsteady and complex flows, has been known for some time. The RANS 
approach assumes that most of the energy is modeled through the turbulence transport equations and is resolved in 
the grid. RANS overpredicts the eddy viscosity, which results in excessive damping of unsteady motion. 
Consequently, the eddy viscosity attains un-physically large values due to unresolved scales, and suppresses most 
temporal and spatial fluctuations in the resolved flow field. Among several methods, the Detached Eddy Simulations 
(DES) [23], the hybrid Large Eddy Simulation (LES) [24-25], the Limited Numerical Scheme (LNS) [26], and the 
Partial Averaged Navier-Stokes (PANS) [27], provide the mechanisms needed to satisfy this requirement. Abdol- 
Hamid and Girimaji [28] explored a new approach to improve the accuracy and robustness of PANS in creating a 
simulation of an unsteady flow field. They accomplished this through the development and implementation of a two- 
stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for 
PANS and other hybrid models. 

This paper describes the history of development of the PAB3D CFD code related to the application of turbulence 
models to jet and nozzle flows, and gives the complete description of turbulence models used in these simulations. 
Several categories of jet flow computation and nozzle configuration [28-64] are briefly described separately in the 
section “PAB3D History of Jet and Nozzle Simulations.” This covers 17 years of the continuous development of 
PAB3D. 


II. PAB3D History of Development 

PAB3D is a structured, multiblock, parallel, implicit, finite-volume solver of the three-dimensional Reynolds 
Averaged Navier-Stokes (RANS) equations, and advanced turbulence models are available in the code. PAB3D is 
widely used for internal and external flow applications by NASA and by the US aerospace industry. PAB3D is one 
of the few codes available today that uses advanced turbulence models in the simulation of complex compressible 
three-dimensional flows. PAB3D developments and applications are documented in over 100 articles. References 
[21, 28, and 29-33] describe some of the major developments in the last 17 years. Some of the capabilities of 
PAB3D are listed below: 

> 3-D RANS/URANS/SMS upwind solver: 

Multiple flow solver (three factor, two factor, and diagonalization) 

Up to third order in space 
Up to second order in time 

> Linear two-equation models with several options: - 

Compressibility correction for high speed flow 
Near wall correction 

Temperature correction for high temperature jet flow 
Up to third order in space 
Up to second order in time 

> Several algebraic Reynolds stress models 

> Tripping option for forced transition from laminar to turbulent flow 

> Multi-scale-type (hybrid) turbulence models: URANS/LES and PANS 

> Steady and unsteady flow boundary conditions 

> Modular multi-block structured with user options selectable for each block: 

Directional grid sequencing 

Turbulence Model options 

Scheme choice including limiter and order 

Viscous terms from simplified to fully coupled 

Number of iterations > 0 (effectively turn block on and off) 

> In-code calculations of integrated forces, moments, and flux quantities 
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> Mixed Roe and van Leer Solver schemes 

> Multi-species and Real Gas Simulation (Frozen Chemistry) 

> Compact memory requirements: average 23 words per grid point 

> GUI capability for input file preparation (UPMS) 

> Conservative patching (PRE) 

> Advanced post processing capability (POST) 

> Advanced jet noise prediction package (JET3D) 

> Dynamic memory allocation 

> Diversified platforms: 

Cray, DEC, HP, IBM, SGI, Linux, Mac and Windows 

> Distributed computing capability 

PAB3D has always been one of the first advanced CFD codes to implement and develop the latest numerical 
schemes, turbulence models, and boundary conditions to enhance its capabilities. Here are some of these examples: 

1- Space Marching Scheme [29], 1989 

2- Multiblock approach [29], 1989 

3- Two-equation turbulence models [12 and 41], 1990 

4- Compressibility corrections [38], 1991 

5- Algebraic Reynolds Stress [34], 1994 

6- High temperature correction [21], 2003 

7- Hybrid methods [28], 2004 


The governing equations of the RANS formulation include the conservation equations for mass, momentum, 
energy, and the equation of state. In the present paper, we use the perfect gas law to represent air properties, and use 
the eddy viscosity and nonlinear concepts to model Reynolds stresses. The mass, momentum, and energy 
conservation equations of the RANS equations can be written in a conservative form as follows: 
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To close the RANS equations, the two-equation (ke) turbulence model is used: 
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C M = 0.09, C el = 1.44 

° k = °k = !-4,o; = (T S = \ andC e2 = C e2 =1.92 

The turbulent stress components are formulated as: 
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- P U i Uj= Ty + Ty L 

The linear contribution to the stress is 

L 2 
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For the purpose of this paper, we will define RANS turbulent viscosity as 
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In the case of linear RANS simulation 


V = v™® 

v t v t 


and C u = 0.09 


Compressibility correction 

PAB3D has several compressibility corrections. In this paper, we will show results based on Sarkar’s 
compressibility correction. 

S = s(l + Ml) 


Where is turbulent Mach number is defined as, M z = . 


3 k 
'2a 


Near wall correction 

PAB3D has several near wall corrections; the most used correction is based on Launder and Sharma defined as: 
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The boundary conditions for k and c at the wall are: 
kwaii = 0 and £ wall ~^ v l 


W 2 
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High temperature correction 

In most applications, we use the Simple Eddy Diffusivity (SED) approach, which is based on the Boussinesq 
viscosity model is used. This approach is used to model all the scalar diffusion terms appearing in the RANS and 
standard k-e equations. For the heat flux term, the SED is written as follows: 


pu t 6 = - 


pv± ar 

(J T dx { 


and g t = 0.9 


For high temperature jet flow, we use the variable developed by Abdol-Hamid et. al. [21] is used as follows: 
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p = 0.09C r 


C T = 1 


where, 
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Algebraic Reynolds Stress 

C /z is 0.09 for the linear model and is a function of vorticity and strain tensors for the nonlinear models. In the 
SZL nonlinear model [35], the turbulent stresses are given by: 
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Stresses in the Girimaji nonlinear model [35] are given by: 
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A compilation of the parameters used in the turbulence models can be found in the Appendix. 

B. Multiscale Approaches 

1 . Two-Stage PANS Approach 

The PANS model [27] was developed to overcome the grid dependency associated with the use of other Hybrid 
Turbulence Models (HTM). In its original form, PANS replaces the two-equation turbulence model by solving for 
the unresolved kinetic energy k u and the dissipation e u . The k u equation is identical to the original k equation. In the 
8 equation (equation 2), the following coefficients are used to change the two-equation model to the Hybrid 
Turbulence Model (HTM), which becomes known as the PANS formulation through the following changes: 

Q 2 = f(cu-Q)+c el 

J £ 

_ f 2 _ f 2 ( 9 ) 

a, - — a, and <x = — a. 


where 
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represent the ratios of the unresolved kinetic energy and dissipation to the total kinetic energy and dissipation, 
respectively. It is only natural to use (/^ , f £ ) to quantify the PANS filter with respect to RANS. Therefore, 

(fk’fs) are used as the resolution control parameters for PANS. Physics of turbulence dictates that: 


(DNS) 0 < f k < f £ < 1 (RANS) 


The original formulation [5] uses constant values for the unresolved kinetic energy parameter (f k ) and unresolved 
dissipation rate parameter (f e ). The users will select values for these parameters and refine the grid until the flow 
solution converges toward a solution target. This could be very time-consuming for resolving complex three- 
dimensional flows. In the present paper, we will discuss an approach to define the unresolved kinetic energy 
parameter. Abdol-Hamid and Girimaji [28] introduced a two- stage approach to estimate the values of the unresolved 
kinetic energy parameter. Here, we will highlight the basic concepts of this approach. Based on a simple 
dimensional analysis, we assume that the turbulent viscosity may be related to the total kinetic energy (k), s, S and A 
as: 


AC 0 ~ £ 

v t » — « AS « A — 


( 10 ) 


which leads to 



3 



Hence, 


fk = c h 





(ii) 


( 12 ) 


C h is a model coefficient, which needs to be calibrated. In the present paper, a value of one will be used to 
evaluate the model. The guidelines, to be followed for the sequential two-stage procedure at this phase of methods 
development, are completely dependent upon flow complexity. Generally speaking, the first stage is a RANS flow 

simulation where the solution serves as the basis to estimate the [f] i ->f £ ) parameters for each grid cell. By using 

this fixed (fi i -,f £ ) distribution over the grid, one would then conduct an unsteady calculation by preferably using a 

high-order scheme and one of the hybrid solver models. An example is demonstrated in the “Unsteady Jet 
Simulation” section near the end of this paper. 


2. Hybrid RANS/LES Approach 

Nichols and Nelson [37] give an example of a hybrid RANS/LES turbulence model. This method was 
implemented in conjunction with Menter’s SST two-equation turbulence model and is termed a multi-scale (MS) 
model. In the present paper, this hybrid model is used with the two-equation model described in equations 2 and 3. 
The turbulent length scale, used in this implementation, is defined as: 

l t = max(6^/ vp NS / Q,/c 3 ' 2 / s ) (13) 
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J1 2 dXj dx , 

_ 1 J 

The sub grid turbulent kinetic energy is defined as: 



The damping function is defined as: 

f d ={ 1 + tanh[2;r(A - 0.5)]} / 2 



( 14 ) 


(15) 


(16) 


X is the unresolved characteristic ratio, and 

A = max(A x , , A z ) 


The eddy viscosity is then calculated from: 


v t =f d vr s Hi-/M ES 

v^=min(v“,0.084A^ 


(17) 

(18) 
(19) 


Note that this hybrid model allows the transition from RANS to LES as a function of the local grid spacing and 
the local turbulent length scale predicted by the RANS model rather than as a function of the grid spacing alone. 
This allows the model to detect whether it can resolve the turbulent scales present on the existing grid before its 
transition over to the LES mode. 


III. PAB3D History of Jet and Nozzle Simulations 

In this section, we will show results for several nozzle configurations simulated using PAB3D code. We will 
target the cases that show the applications of the following major additions to PAB3D: 

1 - Space Marching Scheme 

2- Compressibility correction 

3- Algebraic Reynolds Stress 

4- High temperature correction 

5- Multiscale Approaches (Hybrid methods) 

A. Supersonic Jet Flow 

This section shows examples of using Space Marching Scheme (SMS) and compressibility correction for 
supersonic jet flows. Figure 1 shows a sketch of the jet flow configuration for on-design axisymmetrical jets. Figure 
2 shows computed centerline velocity profiles for a M = 2.22 jet and the experimental data measured by Eggers. 
(65.) The solutions obtained by using different compressibility corrections are compared with experimental data. 
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With no compressibility correction, the potential core length is underpredicted. The location of the end of the 
potential core appears to agree with the centerline velocity profile predicted using the Wilcox model. However, the 
potential core length of the Sarkar solution LJR = 27.15 agrees very well with the data obtained farther downstream. 
The agreement between computational and measured data is much better when compressibility corrections are 
applied, although a small difference exists between the Sarkar model and the Wilcox solutions. Figure 3 shows the 
corresponding results of velocity distributions in the jet cross section at x/R = 25. The importance of compressibility 
correction for supersonic jets is further illustrated here, as the compressibility-corrected computations come very 
close to the measured data, whereas the uncorrected computation underpredicts the centerline velocity by nearly 40 
percent. 

Figure 4 details computed centerline pressure distributions and experimental data for a Mach 2.0 jet at p e lp Q = 
1.445 (Nozzle Pressure Ratio, NPR = 11.3). Three solutions are obtained by using the basic k-z turbulence model 
with no compressibility correction, the Sarkar correction, and the Wilcox correction. The solution without 
compressibility correction shows that the amplitude of pressure oscillations diminishes rapidly downstream and the 
predicted wavelength is much shorter than the experimental data in the downstream region of the jet. The solutions 
obtained with a compressibility-corrected k-z turbulence model show general agreement with measured data. 
Differences between the Sarkar and Wilcox corrections are small. The amplitudes of the computed solutions closely 
follow the test data, but their phase relations with respect to the measured data are somewhat different. At x/R = 40, 
the Wilcox solution leads the measured data by approximately one sixth of one period, whereas the Sarkar solution 
lags behind the measurements by approximately half that amount. All three solutions are very similar near the jet 
exit. However, the amplitude of the first pressure peak near the jet exit is underpredicted by approximately 15 
percent. 

For most jet flows where strong shocks are absent in the computational domain, the space-marching solver in the 
PAB3D code can be used. When the space-marching option is used for jet flow computation as conditions permit, 
the computer time is one-twentieth of the time required for obtaining a time-marching solution with the same flow 
conditions. The accuracy of the solutions obtained by these different solvers is practically indistinguishable. 
Substantial savings in computer time can be realized by using the space-marching method in the PAB3D code if the 
analyses of many cases of jet flow conditions are required for design applications. 

B. Fluidic Injection Nozzle 

The objective of the study was to investigate the performance of multiple injection ports in a two-dimensional, 
convergent-divergent nozzle for fluidic thrust vectoring through experimental and computational means. In 
particular, a two-dimensional, convergent-divergent nozzle with two injection ports was hypothesized to enhance the 
thrust vectoring ability of same nozzle with a single injection port without increasing the total secondary flow 
requirements or incurring significant performance penalties. In addition to verifying this hypothesis, it was the goal 
of the study to validate the computational fluid dynamics code, PAB3D, for multiple injection ports in a two- 
dimensional, convergent-divergent nozzle. The validation of the computational fluid dynamics code would allow 
future investigations to be performed solely by computational means. The test nozzle was designed with 
interchangeable divergent flaps in order to investigate the effects of varying the distance between two injection ports 
in the upper divergent section. A total of five multi-port injection configurations were tested. The divergent flap 
geometry with multi-port injection and a table of important upstream injection port geometry parameters for each 
configuration. This case was chosen to demonstrate the ability of Algebraic Reynolds Stress model in predicting 
complex flow. 

A computational study was conducted on two of the experimental configurations: configuration 1 and 
configuration 5. The configurations were run at NPRs of 4.6, 7, 8.78, and 10, all with a secondary pressure ratio, 
SPR of 0.7. A total of seven boundary (see Figure 5) and initial conditions were used to properly constrain and 
initialize the flow. The initial conditions consisted of a Mach number applied to the inflow duct of the nozzle and to 
the secondary flow plenum to initialize the flow. A wall ’’trip" point, located near the beginning of the inflow duct, 
was used to initialize the turbulent boundary layer of the flow. A stagnation condition was applied to the inflow duct 
of the nozzle and to the plenum of the secondary flow. The stagnation condition was chosen to match experimental 
conditions for total temperature and pressure for the inflow duct and secondary plenum. The static ambient region 
surrounding the nozzle was defined by various conditions. First, the left face was defined by a subsonic inflow 
condition. On the top and bottom faces, a characteristic boundary condition was defined. Finally, a smart boundary 
condition that switched between constant pressure outflow (subsonic) and first order extrapolation (supersonic), 
depending on the local Mach number, was applied to the right face. All solid walls were treated as no-slip adiabatic 
surfaces. 
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Figure 6 shows the predicted flow characteristics using CFD Schlieren flow visualization for configuration 5, 
NPR=4.6, SPR=0.7. Qualitatively, PAB3D predicted the shock patterns very well when compared with the 
experimental flow patterns. Figure 7(a-d) shows the experimental and computational centerline pressures for 
configuration 5, SPR=0.7 (4% of primary mass flow rate). The PAB3D-predicted static pressures along the upper 
and lower nozzle surfaces correlated well with experimental data at all NPRs with a few notable exceptions in shock 
location, strength, and profile. The first shock location on the upper surface at all NPR’s, and the shock location on 
the lower surface at NPR=4.6, was predicted slightly downstream of experimental data. In particular, PAB3D 
predicted the first shock at x/x t =1.3 on the upper surface compared to a shock location of x/x t =1.24 in the 
experimental data at NPR=4.6 (Figure 7(a)). Where x t is the axial location of the throat. Also shown in Figure 7(a), 
PAB3D predicted the shock at x/x t =1.95 on the lower surface compared to a shock location of x/x t =1.89 in the 
experimental data. Second, PAB3D predicted a slightly stronger first shock on the upper surface compared with 
experimental data for all NPR, as indicated by the larger pressure rise in figure 7. On the lower surface, the start of 
the shock is present only at NPR=4.6. A complete comparison of experimental and predicted pressure rise is not 
possible. Finally, PAB3D predicted a steeper pressure profile of the first shock on the upper surface at all NPRs, 
when compared to experimental data. 

Results of the computational study indicate that PAB3D is useful tool for predicting the internal performance 
quantities of a two-dimensional, convergent-divergent nozzle with multiple injection ports. For many subsequent 
secondary flow injection nozzle studies, PAB3D was used with confidence in exploring different configurations for 
multiple injection ports. 

C. High Temperature Jet Flow 

The purpose of this study was to provide insight into the flow physics of pylon-jet interactions and to provide an 
input mean flow field to acoustic prediction methods for installed jet noise. The flow solution needs to be highly 
accurate because it is the basis for discovering important flow phenomena and it is used as input to acoustic 
prediction models. Four configurations were considered; a baseline round nozzle with and without a pylon, and a 
chevron core nozzle with and without a pylon. 

The baseline configuration is a separate flow nozzle with a bypass ratio of five and with an external plug. The 
nozzle and pylon designs are from a nozzle study performed by Boeing in 1996 and represent a generic design. The 
chevrons were designed for the core nozzle using guidelines similar to those used in the NASA Advanced Subsonic 
Transport program. The chevrons were designed to penetrate into the core flow by approximately the thickness of 
the boundary layer. The four configurations analyzed in this study are: 

Configuration 1 : Baseline round core and fan nozzle with no pylon (Figure 8a). 

Configuration 6: Baseline core and fan nozzle with pylon (Figure 8b). 

Configuration 3: Eight-chevron core nozzle and baseline on nozzle with no pylon (Figure 8c). 

Configuration 4F : Eight chevron core nozzle and baseline fan nozzle with pylon and with the tip of a 

chevron aligned with the pylon (Figure 8d) 

For this study, a temperature correction turbulence model, which was tuned on the axisymmetric grid to the 
round nozzle experimental results, was used. Particular attention was paid to matching the axial total temperature at 
x=D = 5 and 10. The model was then frozen for the other cases. Computed total temperature along the axis is plotted 
with experimental data in Figure 9. Comparing the data points at x=D of 5 and 10, configuration 1 is on data, 
configuration 3 is less than 1 D under mixed, and configuration 4F is between 1 and 2 D under mixed. Relative 
errors compared to the predicted round core length are less than 11% for configuration 3, and 11% to 22% for 
configuration 4F. 

The effect of the temperature gradient modified k-s model over the standard k-s model for the round nozzle is to 
increase turbulent mixing and hence decrease the potential core length by 3.4 core diameters, D. The core lengths 
are defined for x=D in which T t = 0.995 T t0 along the axis of the jet; see Figure 9. The action of the modified ks 
model in terms of the increased eddy viscosity closure coefficient is shown in Figure 10, where the ratio of the 
actual to standard value of 0.09 is plotted. In the code, the maximum ratio is set to 5; however, its values in most of 
the flow fields are well below 3. As a result of the undisturbed shear layers, which create strong temperature 
gradients, the round jet configuration exhibits the highest amount of C ^ amplification. CFD solutions of the total 
temperature field are depicted for the round nozzle, Figure 11, with the standard ks model followed by all 
configurations with the modified ks model. 
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D. Unsteady Jet Simulation 

The hybrid RANS/LES and PANS turbulence models are relatively new, and will need to be exercised for a 
wide variety of problems to determine their accuracy before they become an accepted tool for fluid dynamics 
modelers. They seem to offer much for unsteady flow applications, but issues such as grid sensitivity need to be 
addressed. Hopefully, more effort will go into these models in the near future, so that they can mature to a 
confidence level suitable for use in everyday applications. The new capabilities have the potential to improve the 
accuracy and robustness for unsteady flow simulation. This new class of turbulence models is inherently grid size 
dependent, since increasing the grid resolution allows smaller and smaller turbulent scales to be resolved. We have 
introduced and implemented a novel two-stage procedure to efficiently estimate the level of scale resolution possible 
for a given flow on a given grid, for Partial Averaged Navier- Stokes (PANS) and other hybrid models. This is a two- 
stage procedure. In the first stage, a RANS simulation with a Standard Turbulence Model (STM) such as k-sis used 
to produce an estimate of f k over the entire grid domain. In the second stage, we supply f k for the selective 
application of a Hybrid Turbulence Model (HTM) such as the PANS formulation in regions where the grid density 
is sufficient to resolve a portion or all of the large-scale flow structures. In the present implementation, f k is a 
function of length-scale and grid size that represents a characteristic length-scale ratio. We selected the subsonic 
high-temperature jet flows to calibrate and validate the PANS approach. This implementation is a first step towards 
adding a variable resolution turbulence model capability to CFD codes. The PAB3D code can now be used to refine 
the PANS formulation and to conduct validation computations using a variety of simple and complex flow physics 
problems. This approach needs to be calibrated, verified and validated for a wide range of flow problems such as 
different temperature jet, cavity and others. 

The present study used the benchmark experiments performed by Bridges [66] at NASA Glenn Jet Flow 
Faboratory with the flow condition as indicated in Table 1, for core, and free stream. In the following table, T t is 
total temperature, P t is the total pressure, and M is the jet Mach number. 


Table 1: Experimental Subsonic Condition 



T,r°Rj 

P t [psi] 

M 

Core 

1400 

17.68 

0.55 

Free Stream 

540 

14.3 

0.01 


The computational mesh is a full three-dimensional grid with 120 cells in the circumferential direction. The 
computational domain is divided into 48 blocks. The super-fine mesh has a total of 4,000,000 cells. Grid points are 
clustered near the solid surfaces and around the shear layer. Uniform streamwise grid spacing for 2 < x/Dj <12 was 
used. The value of y + for the first cell off the surface varied between 0.2 and 2. During the course of simulating this 
case, single and double precision calculations are used. No significant difference in the results was found. 4, 8 and 
12 sub-iterations for the dual step time accurate approach were used. We also found no significant difference in the 
results. Based on these results, single precision and 4 sub-iterations were used in all of the presented results in this 
paper. 

The computations used the RANS formulation to get time-averaged quantities to calculate the characteristic 
length- scale ratio. This ratio varies in space and is used to produce the unresolved kinetic energy parameter (f k ). 
Figure 12a shows the distribution of this function. This parameter identifies the RANS and PANS regions. The 
RANS regions are defined with the parameter set at a value of one. The PANS regions are the remaining flow 
domains where Rvalues are less than one. The next stage of computation uses this parameter in the execution of the 
PANS formulation. Calibration of the constant C h used the medium grid (1,000,000 cells) to get the velocity profile 
to closely resemble the experimental data for the velocity distribution. It was determined in this process that C h 
should be 1.05 for the present test case. All the calculations performed hereafter used the same value of Ch=1.05. 
Figure 12b shows the two-dimensional snapshot of a slice on the X-Y plane for the velocity (u/UJ) using the PANS 
formulation. The result shows the unsteady behavior of the jet flow as it interacts with the external flow. A similar 
observation was found for a RANS/FES formulation. In the case of RANS prediction, there was not a significant 
difference between the snapshot and the time-averaged flow quantities. This was caused by the fact that RANS over- 
predicts the eddy viscosity, resulting in excessive damping of unsteady motion. Consequently, the eddy viscosity 
attains un-physically large values due to unresolved scales, and suppresses most temporal and spatial fluctuations in 
the resolved flow field. 

Figure 13 shows a comparison of the jet centerline velocity normalized to nozzle exit velocity using RANS, 
RANS/FES and PANS formulations as compared with experimental data. Both RANS/FES and PANS produced 
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better results as when compared with RANS solution. In the case of the RANS results, the resolved kinetic energy, 
k r is zero. Both RANS/LES and PANS provided a mechanism for the RANS equations to resolve the largest scales 
of motion. Figure 14 shows the reduction of the unresolved kinetic energy, k u , as compared with total value. 

IV. Concluding Remarks 

The many steps of development of the PAB3D code for jet flow analysis have demonstrated that Navier-Stokes 
equations with properly implemented two equation k-s turbulence closure models can accurately predict the mixing 
jet exhaust flow from simple or complex three dimension nozzles. The early PAB3D computations had duplicated 
the results of axisymmetric subsonic jet simulations using methods based on the mixing length hypothesis and one 
dimensional momentum theory that were common in the 1970's and into the early 1980’s. Computations for 
supersonic jets had confirmed that code can also predict jet shear mixing in compressible flow when compressibility 
corrections were added to the turbulence model, and that the computed solutions for shock and Mach disk 
containing jets compared very well with classical data in the published literature. By putting jet exhaust flow 
computation on the basis of RANS method and published turbulence models, which are the same as those used for 
aerodynamic computations, it became feasible to analyze internal nozzle performance of complex three dimensional 
nozzles including multi-stream configurations, chevron nozzles, and internal mixers. The applications had further 
extended to the parametric study and design of fluidic nozzle, thrust vectoring, and propulsion airframe integration. 
Another important development was the application of PAB3D code to propulsion airframe aeroacoustics. In 
previous applications, the turbulence models were judged mainly on their ability to ensure an accurate prediction of 
mean flow properties and the performance coefficients of propulsion components. For aeroacoustics applications, 
however, the predicted turbulent kinetic energy and turbulent stresses are applied directly for noise prediction 
purposes. Hence, the standard of accuracy has suddenly elevated to a new level and perhaps to an unknown 
territory. Although early predictions using PAB3D and the JET3D code show great promises, the adequacy of the 
entire analytical basis for the RANS approach and various turbulence model formulations are again open to 
questions and intense scrutiny. In this process, a temperature correction to the turbulence model was implemented 
and tested for hot jet applications, and a broad range of hybrid Navier Stokes computations methods were examined 
for their relevance to modeling the fundamental unsteadiness of jet flow and their application to jet aeroacoustics 
predictions. It is fair to say at this point that many years of interesting and fruitful research in the area are still ahead 
in the future. 
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Appendix Algebraic Stress Models Equation Parameters 

The following functions and variables were used in the algebraic Reynolds stress models used by the SZL 
model: 
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The following functions and variables were used in the algebraic Reynolds stress models used by Girimaji 
model: 
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Figure 1. Typical on-design jet flow configuration and 
terminology. 
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Figure 2. Centerline velocity distribution for supersonic jet using k-z 
turbulence model with different compressibility corrections. 
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Figure 3. Axial velocity component distribution along radial 
direction at x/R = 25 using k-E turbulence model with different 
compressibility corrections. 



Figure 4. Centerline pressure distribution computed with Jones-Launder 
k-E model with different compressibility corrections for pe Ipo = 1.445. 
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Figure 5. Boundary Conditions. 



Figure 6. CFD Schlieren for configuration 5, 
NPR=4.6, SPR=0.7. 
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(b) NPR= 7. 





(c)NPR=8.78. 



Figure 7. Experimental and computational 
centerline pressures for configuration 5, SPR=0. 7. 
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a) Coufig. 1 


b) Coufig. 6 



c) Coufig, 3 d ) Coufig, 4F 


Figure 8. Surface and symmetry plane grid lines. 



Figure 9. Total temperature along jet axis. 
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Figure 10. CFD: Cu factor on symmetry plane. 
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Figure 11. CFD: Total temperature on symmetry plane. 
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a) Unresolved Turbulent Kinetic Energy 



b) Snapshot of Velocity Contour 

Figure 12. Single Nozzle Predictions using PANS Formulation 
Parameter, f k Contours Plane. 
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Figure 13. Normalized Jet Centerline Velocity Results using RANS, 
RANS/LES and PANS Formulations as Compared with 
Experimental Data. 



Figure 14. Normalized Unresolved Turbulent Kinetic Energy, 
k u as Compared with Total Kinetic Energy, k t using PANS. 
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